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scnPtr = img + (S. ?f sM+S . yWinM- 1 ) *pixPerScri S.xOfsM; 
if (levels==8) { 

if (S . segmntMode==3 ) { ./* directed acuity mode 3; pack pixels */ 
PFTsTDTY A /* S.xWinM and S.yWinM known to be mod 8 & mod 2, respectively */ 

for (j-0; j<S.yWinM; j+=2, scnPtr-=2*pixPerScn) { 
pxlPtr = scnPtr; 

for (i=0; i<S.xWinM; i+=2, pxlPtr+=2) { 
bA = (int) *pxlPtr; 
bB = (int) * (pxlPtr+1) ; 
bC = (int) * (pxlPtr-pixPerScn) ; 
bD » (int) * (pxlPtr-pixPerScn+1) ; 

/* bit 0 is set if halftone (&1) */ 
/* bit 1 is set if verticle (&2) */ 

/* code for diagionals is same value in 224 and 28 Iocs */ 
/* diagionals made only when bit d is blacker than a,b,c */ 
vector = 0 ; 

halftone A * / 
halftone B */ 
halftone C */ 
halftone D */ 



0 
□ 

m 
m 
m 



ru 
- 



= 8; 

= 4; 

= 2; 

= 1; 



if (bA&l) vector 
if (bB&l) vector 
if (bC&l) vector 
if (bD&l) vector 
switch (vector) { 

case 0-.lt/* ( ( (bD&224) < (bB&224) ) && ( (bD&224) < (bC&224) ) ) */ (0) 
byt = (bD&227) | ( (bD>>3) &28) ; break; 
} else { 

byt = (bA&2)? (bA&227) | ( (bB>>3) &28) 

: (bA&227) | ( (bC>>3) &28) ; 
byt = grad(byt) ; break; 



case 1 : 



case 2 : 



case 5: 



case 6 : 
case 7: 



byt = (bA&2)? (bA&227) I ( (bB>>3) Sc28) 

: (bA&227) I ( (bC>>3) &28) ; 
byt = grad(byt); break; 

byt = (bB&2)? (bA&224) j ( (bB>>3) Sc28) | (bB&3) 
: (bA&224) jm28 (bA) | (bB&3) ; 

byt = grad(byt); break; 

byt = (bA&2)? (bA&227) I ( (bB>>3) &28) 
: (bA&227) |m28 (bA) ; 

byt = grad(byt) ; break; 

byt = (bC&2)? (bA&224) Im28 (bA) | (bC&3) 

: (bA&224) I ( (bC>>3) &28) | (bC&3) j 

byt = grad(byt); break; 

byt = <bA&2)? (bA&227) Im28 (bA) 

: (bA&227) | ( (bC>>3) &28) ; 

byt = grad(byt); break; 

byt = bA&227 Im28 (bA) ; break; 

byt = bA&227 |m28 (bA) ; 

byt = grad(byt) ; break; 



case 8:if/* ( ( (bD&224) < (bB&224) ) ( (bD&224) < (bCSc224) ) ) */ (0) 
byt = (bD&227) | ( (bD>>3) &28) ; break; 





} 


else 










byt = 


(bD&2)? m224 (bC) 1 ( (bB>>3) &28) 


(bD&3) 






: m224 (bB) 1 ( (bC>>3) &28) 


(bD&3) 




} 


byt = 


grad(byt); break; 




case 


9 : 


byt = 


(bB&2)? m224 (bC) I ( (bB>>3) £c28) 


(bB6c3) 






: m224 (bB) | ( (bC>>3) &28) 


(bB&3) 






byt = 


grad(byt); break; 




case 


10 : 








case 


11 : 


byt = 


m224 (bB) | ( (bB>>3) &28) 


(2) ; 






byt = 


grad(byt) ; break; 




case 


12 : 








case 


13 : 


byt = 


m224 (bC) | ( (bC»3) &28) 
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byt = grad(byt); break; 
case 14 byt = (bD&227) | ( (1 >3)Se28); break; 

case 15: byt = bA; break; 

tempi = (uch)byt; 

putc (tempi, stmPtr) ; /* even is lineart */ 

} else { /* mode must be 0,1 or 2 */ 
/* S.xWinM known to be mod 4 */ 
for (j=0; j<S.yWinM; scnPtr-=pixPerScn) { 

pxlPtr = scnPtr; 

for (i=0; i<S.xWinM; i++, pxlPtr++) { 
tempi = *pxlPtr; 

putc (tempi, stmPtr); /* even is lineart */ 

} else { /* dont do segmentation for quants less than 8 */ 
/* quantize the value, then map into a byte, 0 - 255 */ 
qdblO = (dbl) (l«levels) ; 

Ia qdbll = (qdblO - .000001) / 255.0; 

jp qdbl2 = 255.0 / (qdblO-1.0); 

H /* S.xWinM known to be mod 4 */ 

fit for (j-0; j<S.yWinM; scnPtr-=pixPerScn) { 

m pxlPtr = scnPtr; 

for (i-0; i<S.xWinM; i++, pxlPtr++) { 
fl tempi = *pxlPtr; 

j if ( (templ==255) | | (templ==0) ) { 

putc (tempi, stmPtr); 
L } else { 

L~: qdbl3 = (dbl) tempi; 

qdbl4 = (floor (qdbl3*qdbll) ) *qdbl2; 
^ putc( ( (unsigned char)qdbl4) , stmPtr) ; 

I 

f close (stmPtr) ; 



nt 

rad(byte) 
nt byte; 
return byte; } 

nt 

radl (byte) 
nt byte ; 

/* returns byte with entrys not equal; le a gradient */ 
int bl, b2, b3, b4 , b5 , vector; 

bl = byte&224; b2 = (byte&28)<<3 
if ( ! (bl==b2) ) { 

return byte; 
} else { 

vector = bl>>5; 
switch (vector) { 

case 7: b3=24; break; /* 224 
case 6: b3=28; break; /* 192 
case 5: b3=24; break; /* 160 
case 4: b3=20; break; /* 128 



*/ 
*/ 
*/ 
*/ 
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pxlc = 
pxld = 
else { 

pxlc = pxla; 
pxld = pxlb; 



( (pxlp+pxlr) /2) j 1; 
( (pxlb+pxrs) 12) 1 1; 
/* r is halftone * 



on 
- 



} else { /* p is halftone - generate pqrs as halftone pixels */ 
pxlp = byteP&254; 

pxlq = (byteQ&D? pxlp : byteQ&254; /* q lineart? */ 

pxlr = (byteR&D? pxlp : byteR&254; /* r lineart? */ 

pxls = (byteS&D? pxlp : byteS&254; /* s lineart? */ 

pxla = pxlp; 

pxlb * ( (pxlp+pxlq) /2) ; 

pxlc = ( (pxlp+pxlr) /2) ; 

pxld = ( (pxlp+pxlq+pxlr+pxls) /4) ; 

*(srcPtr+0) = (uch)pxla; 
*(srcPtr+l) = (uch)pxlb; 
*(ofsPtr+0) = (uch)pxlc; 
*(ofsPtr+l) = (uch)pxld,- 

srcPtr += 2; /* advance by one source pixel */ 
ofsPtr += 2; /* advance by one source pixel */ 
nxtPtr += 2; /* advance by one source pixel */ 
IstPtr += 2; /* advance by one source pixel */ 

srcPtr += xwin*2; /* advance by one more source raster */ 

ofsPtr += xwin*2; /* advance by one more source raster */ 

nxtPtr += xwin*2; /* advance by one more source raster */ 

IstPtr += xwin*2; /* advance by one more source raster */ 
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case 3 : b3=8;/ 

case 2: b3=4; 

case 1: b3 = 0; 

case 0: b3=4; 

b4 = byte&3 ; 
b5 = bl|b3|b4; 
return b5; 
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bfeak; /* 96 */ 
break; /* 64 */ 
break; /* 32 */ 
break; /* 0 */ 



it 

.8 (byte) 
it byte; 

if (byte&128) return 28; else return 0; 



it 

224 (byte) 
it byte; 

H if (byte&128) return 224; else return 0; 

5 

m 
m 
H 



©1998 Xerox Corporation 



54 



-cde = = 2. • /* .sasse^le das bytes into .c,z,± pixels * 
srcPtr = srcStrt; 

ofsPtr = srcStrt + xwin*2; /* offset by one raster */ 
nxtPtr = srcStrt + xwin*4; /* offset by two rasters */ 
IstPtr = srcStrt - xwin*4; /* offset by two rasters */ 
for (i=0; i<ywin; { 
for (j-0; j<xwin; { 

byteP = *srcPtr; 
byteQ = * (srcPtr+2) ; 
byteR = *nxtPtr; 
byteS = * (nxtPtr+2) ; 

byteU = *lstPtr; /* U & V have been unpacked last raster */ 

byteV - * (srcPtr-2) ; 

crnrP = * (lstPtr-2) ; 

crnrQ = * (lstPtr+2) ; 

crnrR = * (nxtPtr-2) ; 

crnrS = * (nxtPtr+2) ; 

dglPS = ( (crnrP&l) &&(crnrS&l) ) ; /* both PS corners linearc */ 
dglQR = ( (crnrQ&l) &&(crnrR&l) ) ; /* both QR corners lineart */ 
/* lsb indicates lineart; next lsb indicates horizontal */ 
if (bytePSd) { /* p is lineart - gen pqrs as lineart pixels */ 

tmpp = (byteP&224) ; 

pxlp = tmpp | (tmpp>>3) ; 

if ( ( (byteQ&2) ~ (byteR&2) ) /* one vert, the other horz */ 

&St(byteQ&l) (byteR&l) /* and both sides are lineart * 

&&(byteU&l) &&(byteV&l) /* and four sides are lineart * 

Sc&UdglPS) || (dglQR) ) ) { /* and either ps or qr lineart 

/* if both diags lineart, replicate lowest along diag. */ 

/* if only one diag lineart, replicate in dir of diag. */ 

pxla = pxlp | 1; 

tbot = byteP&28; 

pbot = tbot | (tbot<<3) | 1; 

if < (dglPS) &&(dglQR) ) { /* both diagionals lineart */ 
/* replicate the lowest */ 

if (pxla<pbot) { /* replicate a to d; b or c interpolated 
pxld = pxla; 

if (byteP&2) { /* P horz; bot of p goes to c; b mtrp 
tmpq = <byteQ&224) ; 
pxlq m tmpq | (tmpq>>3) ; 
pxlb - ( (pxlp+pxlq) /2) |1; 

pxlc = pbot ; m 
) else { /* p vert; bot of p goes to b; c intrp */ 
tmpr = (byteR&224) ; 
pxlr - tmpr | (tmpr>>3) ; 
pxlc = ( (pxlp+pxlr) /2) |l; 
pxlb - pbot; 

} else { /* replicate b to c; d interpolated */ 
pxlb = pbot; 
pxlc = pbot; 
tmpq = (byteQ&224) ; 
pxlq - tmpq | (tmpq>>3) ; 
tmpr = (byteR&224) ; 
pxlr = tmpr | (tmpr>>3) ; 
pxld = ( (pxlq+pxlr) /2) |1; 

} else { /* only one diagional PS or QR is lineart */ 
if (dglPS) { /* only diagional PS is lineart */ 
/* replicate along PS direction */ 
pxld = pxla; 

if (byteP&2) { /* p is horz */ 
pxlc = pbot ; 
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p • i = ~mpq; i Cmpq> > 3 ) ; 

pxlb = ( (pxlp+pxiq) /2) | I; 

} else { /* p is vert */ 

pxlb = pbot; 

tmpr = (byteR&224) ; 

pxlr = tmpr | (tmpr>>3) ; 

pxlc = ( (pxlp+pxlr) 12) |1; 

} else { /* only diagional QR- is lineart */ 
/* replicate along QR direction */ 
pxlb = pbot; 
pxlc = pbot; 
tmpq = (byteQ&224) ; 
pxlq = tmpq | (tmpq>>3) ; 
tmpr = (byteR&224) ; 
pxlr = tmpr | (tmpr>>3) ; 
pxld = ( (pxlq+pxlr) /2) |l; 

i=i; /* breakpoint */ 
t } else { 

w pxla = pxlp|l; 

O if (byteP&2) { /* p is horz */ 

fIJ tmpc = byteP&2 8; 

ill pxlc = tmpc | (tmpc<<3) 1 1; 

m if (byteQ&l) { /* q is lineart */ 

SI tmpq = (byteQ&224) ; 

U pxlq = tmpq | (tmpq>>3) ; 

, if (byteQ&2) { /* q is horz */ 

L tmqs = byteQ&28; 

pxqs = tmqs | (tmqs<<3) ; 
} else { /* q is verticle */ 



f.t } else t /» q is vet uj.<_.l<= -/ 

jif if (byteS&l) { /* s is lineart */ 

i; tmps = (byteS&224) ; 

O pxls = tmps | ( tmps>>3 ) ; 

H= pxqs = (pxlq+pxls) /2; 

} else { /* s is halftone */ 

^ pxqs = pxlq; 

pxlb = ( (pxlp+pxlq) /2) I 1; 
pxld- ( (pxlc+pxqs) /2) |1; 
} else { /* q is halftone */ 
pxlb - pxla; 
pxld » pxlc; 

} else { /* p is vert */ 
tmpb * byteP&28; 
pxlb » tmpb | ( tmpb< < 3 ) 1 1 ; 
if (byteR&l) { /* r is lineart */ 
tmpr = (byteR&224) ; 
pxlr = tmpr | (tmpr>>3) ; 
if (byteR&2) { /* r is horz */ 

if (byteS&l) { /* s is lineart */ 
tmps = (byteS&224) ; 
pxls = tmps | (tmps>>3) ; 
pxrs = (pxlr+pxls) /2; 
} else { /* s is halftone */ 
pxrs = pxlr; 

} else { /* r is verticle */ 
tmrs « byteR&2 8; 
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pxlc = ( (pxlp+pxlr) /2) I 1 ; 

pxld = ( (pxlb+pxrs) /2) | 1; 

else { /* r is halftone */ 

pxlc = pxla; 

pxld = pxlb; 



else 
pxlp 
pxlq 
pxlr 
pxls 
pxla 
pxlb 
pxlc 
pxld 



{ /* p is halftone - generate pqrs as halftone pixels 
byteP&254; 

q lineart? */ 
r lineart? */ 
s lineart? */ 



(byteQ&l) ? pxlp 
(byteR&l) ? pxlp 
(byteSSd) ? pxlp 

pxlp; 
( (pxlp+pxlq) /2) ; 
( (pxlp+pxlr) /2) ; 
( (pxlp+pxlq+pxlr+pxl 



byteQ&254 ; 
byteR&254; 
byteS&254; 



) /4) ; 



*(srcPtr+0) = (uch)pxla; 
*(srcPtr+l) = (uch)pxlb; 
*(ofsPtr+0) = (uch)pxlc; 
*(ofsPtr+l) = (uch)pxld; 

srcPtr += 2; /* advance by one source pixel */ 
ofsPtr += 2; /* advance by one source pixel */ 
nxtPtr += 2; /* advance by one source pixel */ 
IstPtr += 2; /* advance by one source pixel */ 

srcPtr += xwin*2; /* advance by one more source raster */ 

ofsPtr += xwin*2; /* advance by one more source raster */ 

nxtPtr += xwin*2; /* advance by one more source raster */ 

IstPtr += xwin*2; /* advance by one more source raster */ 
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